function [quad,R =Romberg(f,a,b，eps)
%[表示被积函数句柄
% a，b表积分区间[a，b]的端点 eps
% eps 表示精度
% quad 是用Romberg加速算法求得的积分值 ead=0
%R为Romberg 表
% err 表示误差的估计 0.10
h=b-a; 0.1;
R(1,1)=h*(feval(f,a)+feval(1i,b))/2; 
M=1;J=0;err=l; 0.
while err>eps
J=J+1;
h=h/2; 
S=0;
for p=]:M
x=a+b*(2*p-1);
S=S= feval(f,x);
end
R(J+1,1)-k(k,1)/2-bS;
M-2M;
for k-1:J
R(J+1,k+1)=R(J+1,k)+(RJ+1;k)-R(J,k))/(4k-1);
end
err=abs(R(J+1,J)-R(J-1，J+1))；end
quad-R(J-1，J+1)；